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We discuss how to simulate simple quantum logic opera- 
tions with a large number of qubits. These simulations are 
needed for experimental testing of scalable solid-state quan- 
tum computers. Quantum logic for remote qubits is simulated 
in a spin chain. Analytical estimates are presented for pos- 
sible correlated errors caused by non-resonant transitions. A 
range of parameters is given in which non-resonant effects can 
be minimized. 

INTRODUCTION 

Several proposals for scalable solidjjstate quantum com- 
puters have been recently publisheairlj. The simulation 
of quantum information processing with a large number 
of qubits will be necessary to test these computers. We 
show in this paper how to implement the Control Not 
(CN) gate between remote qubits in a spin chain using 
a sequence of electromagnetic pulses. The errors gen- 
erated by non-resonant effects are estimated analytically 
and computed numerically. We give the probability of er- 
ror as a function of the Rabi frequency and the number 
of qubits. 

The two-level approximation for an individual transi- 
tion used in this paper for analysis of the system (see also 
Ref. [4]) allows one to simulate the quantum dynamics 
in a system with an enormous number (2 L , where L is 
the number of qubits) of quantum states. This technique 
is based on selection of states generated as a result of the 
resonant or near-resonant transitions, while other transi- 
tions are neglected. 

In Sec. I we give the Hamiltonian of the nuclear spin 
quantum computer. In Sec. II we use the resonance ap- 
proximation to show that each quantum state interacts 
resonantly only with a single state. This allows us to 
decrease from exponential in L to polynomial in L the 
number of differential equations in the system. In Sec. 
Ill we use the solution of the problem in the resonance 
approximation to estimate error probabilities in the im- 
plementation of a CN gate caused by finite detuning from 
exact resonance. 

I. DYNAMICS OF THE SPIN CHAIN 

We consider a chain of identical nuclear spins placed in 
a nonuniform high external magnetic field B, with a uni- 
form gradient. The nuclear magnetic resonance (NMR) 
frequency for the fcth spin is lu^ — j n Bk, where 7„ is the 
nuclear gyromagnetic ratio and is the z-component of 
the magnetic field at the location of the fcth spin. The 
gradient of the magnetic field (in the direction of the 



chain) provides a shift of Uk by the value 8u> between 
the neighboring spins. (For the physical parameters see 
Ref. [4].) 

The Hamiltonian of the spin chain in an external radio- 
frequency (rf) field is, 

L-l L-l 

H„ = — 



k=0 



u k Ii -2JJ2 Iklk+i +V n = H + V n , (1) 
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where J is the Ising interaction constant and Iji. 
the operator of the z-component of spin 1/2. The 
operator, V n , describing the interaction of spins with 
the rf pulses, can be written in the form,tj V n = 

-(fi n /2) Y.kZlVk ex P {-i{ v nt + fn) + I~k exp(w n t+ 
ip n )], where Cl n , v n and (p n are the Rabi frequency, 
the frequency and the phase of the nth pulse; = 
1% ± I\. Each quantum state of the spin chain can be 
described as a superposition of eigenstates of Hq, for ex- 
ample, 1 00 . . . 00), 1 00 . . . 01), and so on, where the state 
| . . . Ofc . . .) corresponds to the direction of the kth nuclear 
spin along the direction of the magnetic field (spin up), 
and the state | . . . 1& . . .) corresponds to the spin being in 
the opposite direction (spin down). The wave function, 
"J, of the spin chain can be written in the interaction 
representation as a linear combination of the individual 
states, \I/(i) = J2 P Q>WIp) ex P( — iE p t), where E p is the 
energy of the state \p) . The contribution of each quan- 
tum state to the wave function, ^(t), is given by the 
coefficient C p (t). The quantity |C p (i)| 2 is the probability 
of finding the spin chain in the state \p) at time t. 

The Schrodinger equation for the coefficients C p (t) 
has the form (we put tp n = 0), iC p (t) = 

Y? m =o V pm ex-p[i(E p - E m )t + ir pm u n t]C m (t), where 
r P m — Tl for E p > E m and E p < E m , respectively. 
Vp l m = —Q n /2 for the states \p) and \m) connected by a 
single-spin transition. V™ — for all other states. 

II. THE APPROXIMATE SOLUTION 

The Schrodinger equation can be numerically inte- 
grated only for a spin chain with small enough number of 
spins since the number of states increases exponentially 
with L. The problem can be simplified when the Rabi fre- 
quency, Q n , is much less than the difference, 5uj, between 
the NMR frequencies of the neighboring spins, il n <C Sw. 
Suppose that a pulse is resonant with the fcth spin in the 
chain, v n ~ u>k- Then one can admit that this pulse af- 
fects in some approximation only the kth. spin in the chain 
and does not interact with other spins. In this case, only 
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the term with the energy E m i ~ E p + v n effectively con- 
tributes to the right-hand side of the Schrodinger equa- 
tion, so that the states \m') and \p) are connected by a 
single-spin flip of the fcth spin. In this case, we reduce 
the Schrodinger equation to the set of two differential 
equations, 

iC m (t) = — ^ exp[i(E p - E m - u n t]C m (t), 

%C p (t) = exp[-i(E p -E m - u n t]C p (t), (2) 

where E p > E m , and \p) and |m) are two stationary 
states which are connected by a single-spin transition of 
the fcth spin, v n ~ ujk- 

The solution of Eq. (Q), for initial conditions, 
Cm(in-i) = 1) C p (t„_i) = 0, after the action of the 
nth pulse has the form, 



Cm i^n ) 



cos (A„t„/2) + i^- sin (A„t„/2) 



Cp(t n ) =z^sin(A„T„/2)e( l *"- lA " +i;; ^), (3) 

where t n _i and t n are the beginning and the and of the 
nth pulse, r„ is the duration of the nth pulse, A„ = 
E p — E m —v n (we do not indicate the dependence of A„ on 
the indices p and m), A„ = yj f2| + A„ is the precession 
frequency in the frame rotating with the frequency v n . 
Suppose that at t n -i the system was in the state \m). 
In the case of the exact resonance (A„ = 0) and for 
QnT n = 7r (a 7r pulse) Eqs. (|^) describe the resonant 
transition from the state \m) to the state \p). If the 
frequency of the next pulse also satisfies the resonance 
condition, A n +i = Ei — E p — v n +i = 0, then the system 
with probability equal to unity will transform to the state 
|Z), and so on. 

III. ERRORS IN CREATION OF ENTANGLED 
STATE FOR REMOTE QUBITS 

The advantages given by the two-level approximation 
described above are obvious: (a) we solve the problem by 
using the discrete map (|^) instead of integrating differ- 
ential equations; (b) we can consider only the states with 
large enough probabilities and neglect all other states in 
a controlled way; (c) we can estimate errors caused by 
non-resonant effects and minimize them by choosing op- 
timal parameters. 

In this paper we will estimate the error in implementa- 
tion of a unitary operation for remote qubits (which is a 
particular case of the well-known Control-Not (CN) gate) 
using a sequence of 7r-pulses. This CN gate is defined as 
the unitary operator, U, with the following properties: 

f/|0 i _ 1 i _ 2 ...0i0 ) =e^ 1 |0 L -i0 L -2...0i0o), (4a) 



C/|1l-iO L -2...0i0 ) 



•Oilo), (4b) 



is in the state |1). The operator, U, can be used to create 
an entangled state for remote qubits: J7(q;|Ol-i---Oo) + 
/?|l L _i...0 o )) =ae^0L-i...0 o )+/3e^|l L _i...l o ). To ac- 
complish the (4b) operation in the system described by 
the Hamiltonian (1), we choose a sequence of 7r-pulscs 
with resonant frequencies (for which all A„ = 0). In 
this case, the state 1 10 . . . 00) transforms to the state 
1 10 . . . 01), with probability equal to unity, by the follow- 
ing scheme: |1000...0) -> |1100...0) -> |1110...0) -> 
|1010... 0) -> 1 1011... 0) -> |1001...0) ->...-> 
1 100 ... 11) -> 1 100 ... 01). The sequence of pulses which 
realizes this protocol has the following form: v\ = UJ2L-2, 

V2 = U 2 L-3, v 3 = LU 2 L-2 ~ 2 J, ^2L-4 = ^6, 

V2L-A = — J, V2L-Z = wi- If we apply the same pro- 
tocol to the state 1 00 . . . 00), then with some probability 
the system will remain in this state because all transi- 
tions are non-resonant with the detuning A„ = ± J, ±2 J. 
Since A„ ^ these non-resonant transitions have the 
probabilities (see Eq. 



(^n/A„) 2 sin 2 (A„T„/2). 



(5) 



In this paper we study the errors generated during 
the operation ( flaj) under the condition that the tran- 
sitions in Eq. ( [lb]) are resonant. The operator U in 
Eqs. ( |4a| ) and ( ffb| ) can be written as a product, U = 
U2L-3U2L-4 ■ ■ ■ UiUx, of operators of individual pulses 
Um where n = 1, , 2, . . . , 2L — 3. We take the Rabi fre- 
quency to be the same, fl n = Q, for all pulses. Then 
the error, e n , generated by the nth pulse is defined only 
by the detuning, A„. The values of detuning, |A„|, are 
the same for all pulses, |A„| = J, except for the third 
pulse, where IA3I = 2 J. We denote e„ = e for n 7^ 3 and 

£n=3 = £ ■ 

Each unwanted state generated by the nth pulse can 
produce other unwanted states under the action of the 
other (n + m) pulses. This causes the generation of a 
hierarchy of different states. Instead of one ground state, 
which should be the result of (4a), we generate a series 
of different unwanted states. 




where ipi and <p2 are known phases. The target qubit (0) 
should change its state only when the control qubit (L—l) 



FIG. 1. Dependence of e (solid line) and e' (dotted line) on 

n. j = 1. 

Since the dynamics is generated by the discrete map 
(^), probability of error can be estimated analytically. 
Suppose that we can measure only the states with prob- 
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ability P > Pq. where Pq can be denned by the exper- 
imental conditions. Then, if e, e' < £1 = Po, we shall 
not find in the system any unwanted states and the CN 
gate is realized with probability equal to unity. In Fig. 1, 
we plot e (solid line) and e' (dotted line) as functions 
of Q,. We suppose that E\ = Pq = 10~ 6 . One can see 
from Fig. 1 that e(Sl) and e'(fi) are less than e\ only in 
very narrow regions of the values CtLf2 (in the vicinity of 
flk = |A„|/V4fc 2 - 1, k = 1,2, ..El). Hence, high pre- 
cision is required to implement the CN gate, (|4a|), (|4q), 
without error. 

Let us consider now in detail the case of arbitrary 
e: < e < 1, and explain how the errors can be es- 
timated analytically in the operation (|ia|). The first 
pulse U\ generates one unwanted state C/i |0000 . . .0) — > 
1 0100 ... 0) with the probability e. The probability to re- 
main in the ground state is 1 — e. The second pulse, 
U2, generates an additional unwanted state from the 
ground state, t/ 2 |0000 . . . 0) -> |0010...0), with prob- 
ability e(l — e), and transforms the state |0100...0) 
to the state 1 01 10 . . . 0) without generating additional 
states, since the second pulse is resonant to this transi- 
tion. The probability of remaining in the ground state is 
now (1 — e) 2 . The total number of states after the action 
of two pulses is three. The probabilities of all generated 
states after the action of the U3 pulse are: 



0000. 


..0): 


(l-e) 2 (l-e 


0100. 


.0): 


e'{l~e)\ 


0010. 


.0): 


e(l-e + e 2 ), 


0110. 


.0): 


e(l-e 2 ). 



(Note, that above expressions can be considered only as 
estimates because we estimated probabilities instead of 
complex amplitudes.) One can continue to calculate the 
states generated by subsequent pulses and estimate their 
probabilities in a similar way. However, for e, e' -C 1, one 
can omit the contribution from the higher order terms in 
e, and neglect the states with the probabilities less than 
Po- 

Suppose, that e' « e and £\ < e < £2 = y/Po- In this 
case, each pulse generates only one unwanted state from 
the ground state. Then, only 

N 1 (L)=2L-Z (6) 

unwanted states will be generated with probabilities 
(with the accuracy up to e 2 ) e, e(l — e), e(l — 2e), . . . , 
e(l — (2L — 4)e). After the action of 2L — 3 pulses there 
are L — 1 unwanted states of the form, 

|011 ... 11), |0011 ... 11), .-, ]00 . . . 011), |00...01), (7) 

with the target spin in the state 1 1) , and L — 2 unwanted 
states of the form, 

I011...110), |0011... 110), |00...010), (8) 

with the target spin in the state |0). As one can see from 
Eqs. (0) and (||), already in the first order of e many 



states have correlated error with many qubits in excited 
states. The estimate of the total probability of unwanted 
states is 

2L-4 , . 

P 1 (£)= S J2(l-ne)^£(l--£ , (9) 

71 = ^ ' 

where £ = (2L — 3)e. 

The probability of unwanted effect of finding a target 
qubit in the state |1) while the control qubit is in the 
state |0) (i.e. the probability of the process U\0 ... 0) — > 
|0...1))is, 

L-3 

Vx{e) = £ + £ ^(1 - ( 2n + ~ r (! - r ), ( 10 ) 

n=0 



where T = (L — 2)s. 
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FIG. 2. (a): The total probability, Pi, of unwanted states 
in the execution of the CN logic gate as a function of the 
number of qubits, L, in the spin chain. Solid line: analytical 
estimate ([]), filled circles connected by dashed line: numerical 
results, (b) Probability, Pi, of observing the target spin in 
the state |1), while the control qubit is in the state JO), after 
the CN operation; solid line: analytical estimate (tOI), filled 
circles connected by dashed line: numerical results. J = 1, 
n = 0.0906. 

In Figs. 2 (a), (b) we plot the probabilities, Pi and Vi, 
as functions of the number of spins, L, in the chain, for 
fl = 0.0906, so that e x < e(fi) = 4.78 x 10~ 5 < e 2 = 1CT 3 
(e'(fi) = 3.23 x 10~ 5 ). The number of unwanted states 
in our calculations was exactly equal to N\ in Eq. (|^). 
From Figs. 2 (a) and (b) one can see that the analytical 
estimates agree with the results of numerical calculations. 

In a similar way one can calculate the probabilities of 
unwanted states for the case £2 < £ < £3 = (Pq) 1 / 3 . In 
this situation one should take into account states with 
the probabilities up to e 2 . Each of these states is gener- 
ated as a result of two successive non-resonant processes. 
However, since the ground state generates only unwanted 
states with the probabilities of order e, the total proba- 
bility of unwanted states will be given again by Eq. (||). 
But the number of generated states increases and at the 
condition, £2 < £ < £3, it is proportional to L 2 . 

In Fig. 3 we compare the numerical results with the 
analytical estimates for the total probability of unwanted 
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states given by Eq. ^ when the value e is relatively 
large. From Fig. 3 one can see that formula (0) correctly 
describes the behavior of the system for large value of e. 
When L becomes large enough, one should also include 
into consideration in Eq. (^) the terms of order (Ls) 3 . 




30 50 70 90 

Number of spins 

FIG. 3. The total probability of unwanted states for the 
case of relatively large e, when ei < e < £2. Solid line: 
analytical estimate (190 , filled circles connected by dashed line: 
numerical results. J = 1, Q = 0.20844, e(fl) = 2.98 x 10~ 3 , 
e'(fi) = 2.41 x 10" 3 . 

The probabilities of unwanted states are shown in 
Figs. 4 (a) and (b) in two different scales. Two types 
of states are clearly seen. The states in Fig. 4 (a) are 
generated as a result of one non-resonant transition, so 
the probabilities of these states are proportional to e. On 
the other hand, the states in Fig. 4 (b) are generated by 
two successive non-resonant transitions, and they have 
the probabilities of the order e 2 . 



(a) 



(b) 
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a large number of qubits. In the limit when the Rabi 
frequency fl is much less than the difference between the 
NMR frequencies of neighboring spins, it is shown that 
the probability of error, P, is mainly defined by the small 
parameter e defined by (5) and by the number of qubits 
in the spin chain. It is demonstrated that at definite 
values of Q, the value of e is small and there is no corre- 
lated error in the implementation of the CN gate. The 
total probability of unwanted states at other values of e 
is estimated analytically and computed numerically. 

Since the Hamiltonian of the system only allows transi- 
tions between the states connected by a single spin trans- 
formation, and because the electromagnetic field inter- 
acts resonantly only with one spin in the spin chain, the 
probability of error cannot accumulate in a definite sin- 
gle unwanted state. The probability of each unwanted 
state is always less than the value of the small parame- 
ter of the problem e(fl). The considered approach can 
be useful for experimental testing of scalable solid-state 
quantum computers. 
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FIG. 4. The probabilities of generated states in two differ- 
ent scales. L — 70. All other parameters are the same as in 
Fig. 3. 



IV. CONCLUSION 

We have demonstrated how one can simulate a simple 
logic operation for a large number of qubits. We analyzed 
the probability of correlated errors in the implementation 
of a CN gate between remote qubits in a spin chain with 
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